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^ ABSTRACT 

IL" We measure the real-space galaxy power spectrum on large scales at redshifts 0.5 

.5^ to 1.2 using optical colour-selected samples from the CFHT Legacy Survey. With 

the redshift distributions measured with a preliminary ^14000 spectroscopic redshifts 
from the VIMOS Public Extragalactic Redshift Survey (VIPERS), we deproject the 
angular distribution and directly estimate the three-dimensional power spectrum. We 
use a maximum likelihood estimator that is optimal for a Gaussian random field 
giving well-defined window functions and error estimates. This measurement presents 
an initial look at the large-scale structure field probed by the VIPERS survey. We 
measure the galaxy bias of the VIPERS-like sample to be bg — 1.38 ± 0.05 (cs = 0.8) 
on scales k < 0.2 h Mpc~^ averaged over 0.5 < z < 1.2. We further investigate three 
photometric redshift slices, and marginalising over the bias factors while keeping other 
ACDM parameters fixed, we find the matter density fi™ — 0.30 ± 0.06. 

Key words: cosmological parameters, observations, large-scale structure of Universe, 
methods: statistical 
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1 INTRODUCTION 

The shape of the galaxy clustering power spectrum encodes 
the dynamical history of the Universe under the influence 
of baryons, dark matter and dark energy. On large scales 
the assumption of Gaussianity can be made, and the statis- 
tic summarises all of the cosmological information that is 
available. 

Measurements of the power spectrum at z ~ have 
led to fundamental tests of the ACDM model (Tegmark 
et al. 2004; Efstathiou et al. 2002). The angular distribu- 
tion of galaxies on the sky, although less sensitive than 
the full three-dimensional view, has played an important 
role as well. Indeed, strong tests of the CDM model were 
made with the two-dimensional correlation function from 
the APM galaxy survey (Maddox et al. 1990). Photometric 
surveys have the capability to probe significantly larger vol- 
umes at higher sampling rates than targeted spectroscopic 
surveys. The advantages have become clear with the ad- 
vancement of photometric redshift estimation methods. The 
loss in three-dimensional precision can be compensated for 
with increased statistics leading to strong cosmological con- 
straints that are comparable to the results from spectro- 
scopic surveys. 

Additionally, the projected density field is only weakly 
sensitive to redshift-space distortions, thus it provides a 
means to infer the real-space power spectrum directly. The 
dependence on peculiar velocities becomes important for 
narrow redshift slices and can be turned into a useful mea- 
sure of the growth rate (Ross et al. 2011a). Measurements 
of the baryon acoustic feature and redshift-space distortions 
have now been made on photometric samples taken from the 
Sloan Digital Sky Survey (Padmanabhan et al. 2007; Blake 
et al. 2007; Thomas, Abdalla & Lahav 2011). 

In this analysis, we present a new measurement of the 
real-space galaxy power spectrum using a photometric cat- 
alogue of galaxies at 0.5 < z < 1.2 from the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS) Wide survey. 
The survey consists of four fields covering a total area of 
133 deg^ The extent of the largest field, Wl, is ~ 10° or 
200 h~^ Mpc at z=0.7, giving a maximum scale we may 
probe of k-min ~ 0.05 h Mpc~^. The dataset has been used 
for previous cosmological analyses, in particular for weak 
lensing (Fu et al. 2008; Kilbinger et al. 2009; Tereno et al. 
2009; Shan et al. 2011) and galaxy correlation function mea- 
surements (Coupon et al. 2011). 

A key ingredient needed to interpret the projected den- 
sity field and constrain the three-dimensional power spec- 
trum is the redshift distribution of the galaxy sample. For 
this we use spectroscopy from the VIMOS Public Extra- 
galactic Survey^ (VIPERS) (Guzzo et al. 2011). VIPERS is 
an ongoing spectroscopic program to target 10^ galaxies in 
the redshift range 0.5-1.2 in a total area of 24deg^ in the 
CFHTLS Wl and W4 fields. The accuracy of the spectro- 
scopic measurements from VIPERS provide an unbiased es- 
timate of the redshift distribution. With this knowledge, we 
are confident that we can deproject the angular clustering 
signal and constrain the three dimensional power spectrum. 
The primary advantage of studying the deprojected power 
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Figure 1. Galaxy count maps in the CFHTLS- Wide fields with 
the VIPERS-like colour selection. We use Healpix cells with size 
7'. Gaps in the survey coverage are left as blank pixels. The grid 
overlay has spacing of 1°. 



spectrum Pk is its closeness to theory. The shape of the an- 
gular power spectrum is complicated by its dependence on 
survey properties, its depth and geometry. Furthermore, in 
projection, scales are mixed. Ideally we would like to sep- 
arate the power on large scales in the linear regime from 
power on small scales that is infiuenced by complex astro- 
physical processes. 

How to derive the three-dimensional power spectrum 
from the two-dimensional density field is a problem of decon- 
volution. A good inversion method should be stable against 
noise in the data. Preliminary work done by Baugh & Efs- 
tathiou (1993, 1994) and Gaztanaga & Baugh (1998) used 
the Lucy deconvolution method that is known to be robust. 
To further derive cosmological constraints, we must be able 
to estimate the covariance of the deprojection. Methods of 
propagating the error from the angular correlation function 
to the three-dimensional power spectrum were developed by 
Dodelson & Gaztaiiaga (2000) who perform the inversion 
with a prior on the smoothness of the power spectrum and 
compute a covariance matrix of the estimate. Further work 
by Eisenstein & Zaldarriaga (2001), Dodelson et al. (2002) 
and Mailer et al. (2005) made use of the singular-value de- 
composition technique to remove modes that destabilise the 
inversion. 

Importantly, the deprojection method should produce 
well-defined window functions that describe the mode mix- 
ing. The aim is to separate the small and large scales that are 
mixed in projection, and the residual leakage should be un- 
derstood. A similar problem was solved with the maximum 
likelihood methods developed for the cosmic microwave 
background angular power spectrum (Tegmark 1997; Bond, 
Jaffe & Knox 1998) and then later applied to galaxy surveys 
(Huterer, Knox & Nichol 2001; Tegmark et al. 2002). Ap- 
plications to the deprojection of the power spectrum were 
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presented by Efstathiou & Moody (2001) and Szalay et al. 
(2003). 

In this work, we adopt the maximum-hkehhood tech- 
nique to construct an estimator for the power spectrum. The 
result is optimal under the assumption that the density is 
represented by a Gaussian random field. This is a reasonable 
assumption for the galaxy distribution on large scales. More- 
over, the estimator also simultaneously gives the covariance 
of the estimate as well as the window functions. For small 
surveys, where the window functions must be handled care- 
fully, the approach is especially useful. We pay close atten- 
tion to the window functions for the results presented here. 
Maximum likelihood estimates are computationally expen- 
sive. However, because of the relatively small field sizes we 
consider, we can perform all computations on a consumer 
level four-core desktop computer. 

In this article, we first introduce the CFHTLS and 
VIPERS datasets used. In Sections 3 and 4 we review the 
angular power spectrum formalism and the maximum like- 
lihood deprojection using a quadratic estimator. We then 
apply the method to Gaussian simulations and investigate 
potential biases due to uncertainty in the redshift distri- 
bution and the fiducial cosmology. Lastly, we measure the 
power spectrum with CFHTLS data and constrain the linear 
galaxy bias and matter density. 

We report magnitudes using the AB magnitude conven- 
tion in the CFHT u* g'r'i' z photometric system. We assume 
a flat ACDM cosmology with J/o = 70.4 km s"^ Mpc"\ 

= 0.272, = 0.0456, = 0.963 and ag = 0.8 (Larson 
et al. 2011). 



2 DATA 

2.1 Photometric selection 

The Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) Wide includes four fields labelled Wl, W2, W3 
and W4. The total area is 133deg^ imaged with five-band 
photometry ugriz to a depth of i = 24.5. We construct 
colour-selected galaxy samples from these fields to match 
the spectroscopic target selection used by the VIMOS Pub- 
lic Extragalactic Survey (VIPERS). 

VIPERS is a spectroscopic program to measure the red- 
shifts of galaxies over an area of 24deg^ in the CFHTLS 
Wl and W4 fields. Galaxies are targeted from CFHTLS- 
Wide photometry to a flux limit of iAB = 22.5 with colour 
criteria to produce a sample at z > 0.5 having few low- 
redshift interlopers. The selection is done in the u — g, 
r ~ i colour plane with the following limits on extinction- 
corrected magnitudes: (1) r — i ^ 0.7 AND u — g ^ 1.4 OR 
(2) r - j > 0.5{u - g) AND u-g <1A (Guzzo et al. 2011). 
We replicate these colour criteria on the full CFHTLS- Wide 
photometry T0006 release (Goranova et al. 2009). Hereafter, 
we refer to this selection as the VIPERS-like sample. 

Every source has an estimated photometric redshift and 
star-galaxy classification from the T0006 photometric red- 
shift catalogue. The star-galaxy classification accounts for 
both the source profile and fits to stellar spectral tem- 
plates (Coupon et al. 2009). We apply the same criteria 
as used for the VIPERS target selection and exclude all 
sources photometrically classified as stars (7% of sources). 
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Figure 2. Redshift distributions of the spectroscopic samples 
used, (a) The distribution of the full VIPERS sample is shown. 
Overplotted are the distributions from the Wl and W4 fields 
individually, (b) The spectroscopic redshift distributions of the 
three photometric redshift subsamples are plotted. 



From the remaining sample identified as galaxies, we remove 
sources that fail to fit galaxy spectral templates with re- 
duced > 100. This cut removes sources with incomplete 
or spurious photometry amounting to ~0.15% of sources 
which are typically near to the edges of masked regions, 
bright stars or field borders. 

We also use the photometric redshifts to select sub- 
samples of narrower slices in redshift, labelled S6 (0.5 < 
zphot < 0.6), S7 (0.6 < zphot < 0.8) and S8 (0.8 < Zphot < 
1.), see Table 1. As with the VIPERS-like sample, these are 
also limited to ias = 22.5. We have confirmed that these 
photometric redshift selections at 2; > 0.5 also meet the 
VIPERS selection criteria. Thus, the VIPERS spectroscopy 
can be used to calibrate the redshift distributions of these 
samples without introducing a bias. 

The VIPERS spectroscopic targets were selected from 
the CFHTLS T0005 catalogue after field-to-field colour cor- 
rections were applied by the VIPERS team. These colour 
corrections are no longer necessary in the T0006 update, 
and it has been shown that the selections from the two cat- 
alogue versions match well (Guzzo et al. 2011). We note that 
a limited area in the CFHTLS has been observed with a re- 
placement i-band filter called y. The photometric redshifts 
were computed with the appropriate filter transmission func- 
tion. To construct our i < 22.5 limited samples, we take a 
reasonable approach and do not distinguish between the two 
bands. 

The CFHTLS catalogues also include corrections for 
Galactic extinction. We do not consider residual systematic 
effects of reddening here because all four fields have rela- 
tively low and uniform extinction at the level of E{B — V) — 
0.06 in W4 and < 0.02 in the other fields. The i = 22.5 limit 
is 2 magnitudes brighter than the detection limit of the sur- 
vey, thus we do not expect the selection to be affected by 
the extinction correction. 
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2.2 Density maps 

Wc construct the density maps for the photometrically- 
selected samples by counting galaxies in cells defined by 
the Healpix scheme with a resolution of 7' (ugide = 512) 
(Gorski et al. 2005). The maps for the VIPERS-Iikc colour 
selection are shown in Fig. 1. We use a survey mask pro- 
vided by CFHTLS and exclude sources that fall within the 
halos of bright stars. For cells that fall on a mask boundary, 
we measure the fractional coverage using a uniformly spaced 
grid of 16 X 16 test points within the cell. Wc use Manglc2 
to test if points fall inside the mask (Swanson et al. 2008). 
This provides us with a weight map Wi = 1/fi where fi is 
the fractional sampling for pixel i. Cells that have less than 
50% inclusion in the survey are removed from the map. The 
areas of the four fields Wl, W2, W3 and W4 are 57.7, 18.6, 
36.8, 19.6 deg^. After putting galaxies in the Healpix cells, 
the number of pixels in the four maps are: 4787, 1592, 3045 
and 1651. 

With Tii galaxies counted in cell i, the over-density is 
computed with 5i = riiWi/n — 1. The mean density in a cell, 
n, is computed from all four fields as n = '^^w^n^/'^^Wi. 
The variance of 5i, assuming Poisson statistics, is erf = 
wf/n. 

The clustering of foreground stars can be a significant 
source of systematic error on large angular scales (Ross et al. 
2011a). For the CFHTLS, we can estimate the stellar con- 
tamination rate independently in each of the four fields and 
apply local corrections to the galaxy density. Wc measure 
the contamination rate directly in the Wl and W4 fields by 
counting the number of targets spectroscopically classified 
as stars in the VIPERS sample. We then extrapolate these 
rates to the W2 and W3 fields by computing the fraction of 
sources photometrically classified as stars and then scaling. 

The total counts in a cell broken down into stars and 
galaxies is given by Nobaerved = Ni,+Ngaiaxy, and the stellar 
contamination fraction is /* = Ni, /Nobaerved- The values we 
derive are listed in Table 2. We apply the correction in the 
following way; Si^corr = 5i/(l-/*) and af ^orr = ^i/i'^-f*)^ 
(Huterer, Knox & Nichol 2001). The effect on the amplitude 
of the power spectrum is ~ 5%. 

A fraction of galaxies will also be misclassified as stars 
and removed from the sample. However, as long as the sam- 
ple is representative of the full population, the power spec- 
trum measurement will not be biased. This may not be the 
case in reality since misclassified galaxies may preferentially 
represent a subclass with a different power spectrum ampli- 
tude. We do not investigate this correction here. 



2.3 Redshift distribution 

Wc use the VIPERS spectroscopic redshift catalogue (inter- 
nal release, version 1.1) to calibrate the redshift distribution 
of the photometric samples, see Table 1. In total, we use 
13191 galaxies from VIPERS including 6516 from the Wl 
field and 6675 from the W4 field. All targets that meet the 
photometric selection criteria with secure redshifts arc used. 
We select based on the quality flag zf lag. Reliable redshifts 
have zflag modulo 10 ^ 2, and we take zflag in {2. .9} 
(galaxy type), {12.. 19} (AGN type) and {22..29} (serendip- 
itous detections). The flag also has a fractional part indicat- 
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Table 2. Star contamination fractions 



Sample 


Wl 


W2 


W3 


W4 


SV 


0.019 


0.056 


0.020 


0.044 


S6 


0.013 


0.030 


0.010 


0.017 


S7 


0.015 


0.042 


0.015 


0.032 


S8 


0.017 


0.055 


0.019 


0.048 



ing agreement with the photometric redshift on a scale from 
1 to 5, where 5 indicates good agreement (within la). 

We estimate the redshift distribution from the his- 
togram of spectroscopic redshifts with a bin size of Az = 
0.05, see Fig. 2. We use the histograms directly in the anal- 
ysis with linear interpolation between bin centres. The red- 
shift distribution of the Wl and W4 samples are remarkably 
similar despite that these fields are well separated on the sky. 
We use the distributions from the two fields to test the im- 
pact of cosmic variance on our results. As we will conclude in 
Section 5.1, small perturbations to the redshift distribution 
do not strongly impact the results. 

The selection function of the VIPERS survey is not uni- 
form with target apparent flux. There are two sampling rates 
that we consider: first, the fraction of potential targets that 
are selected for observation, and second, the fraction of ob- 
served targets that give a successful redshift measurement. 
The first distribution is nearly uniform; the VIMOS spec- 
trograph can place slits on ~40% of the potential targets 
and this fraction is found to be independent of the magni- 
tude of the source. However, we do find that the fraction of 
targets that have a measured redshift with qualities meeting 
our criteria drops from 100% at i = 19 to 50% at i = 22.5, 
the flux limit of the survey. This trend may be corrected for 
by weighting the contribution of each galaxy in the redshift 
distribution by the inverse of the sampling rate. However, 
we find that the correction has a negligible effect on the dis- 
tribution. Weighting the galaxies shifts the mean redshifts 
of the samples by less than Az = 0.01. We also confirm that 
lowering the quality threshold of the spectroscopic sample to 
zflag >= 1.5, which adds 11% additional sources, does not 
significantly alter the distribution. In Section 5.1 we con- 
sider the effect of shifting the mean redshift by Az = 0.05 
to provide an overly-conservative check on the effect of un- 
certainties in the redshift distribution. 



3 ANGULAR POWER SPECTRUM 

From galaxy counts in an imago we may infer the 
projected over-density of galaxies on the sky, S{n) = 
/q°° <53D(n, r)0(r)r^dr. Typically this is an integration 
through a broad slice in redshift defined by a photomet- 
ric galaocy selection function or simply by the limiting flux 
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Figure 3. The correlation matrix for Pj; estimated from the 
VIPERS-Uke sample. Elements of the matrix are labelled with 
the per-cent correlation. Although the window functions overlap 
significantly, the bins are nearly statistically independent by con- 
struction. 



of the survey. We expand the density field in spherical har- 
monics and express the power in mode I by the spectrum 
C. 

We may write the angular power spectrum as a pro- 
jection of the three-dimensional power spectrum, Pu — 
(|^3i3,*:|^). On large scales, the power spectrum evolves with 
the linear growth factor, Di(z). We scale the power spec- 
trum taken at the median redshift of the sample, z, as 
Pfe(z) = [Di{z)/Di{z)f Pk(z). This gives, 



Ci = 



{r)Di{z)/Di{z)ji{krydr 



dk 



(1) 



where r is comoving distance. In the small angle approxi- 
mation, the spherical Bessel function can be approximated 
as ji(x) = y^jiqri + | ^ ^)) expressed with a Dirac delta 
function, and we find Limber's equation: 



Ci = J gi{k)Pk 



dk 
~k' 



(2) 



9i{k) = 



[r'^ir)D^{z)/D^{z)Y 



at r = 



1 + 1/2 
k ■ 



(3) 



The projection kernel, giik), is given by, 
1 

1 + 1/2 

A correction may be added to account for redshift space dis- 
tortions although it is sizable only on large scales at / < 50 
that we are not sensitive to here (Thomas, Abdalla & Lahav 
2011; Ross et al. 2011b). We may now approach the depro- 
jection problem as a deconvolution of Limber's equation. 



4 POWER SPECTRUM ESTIMATOR 

On large scales, the galaxy density may be described by a 
Gaussian random field and the distribution is fully charac- 
terised by its variance. With this assumption, the likelihood 
function of the observed overdensities on the sky may be 
written explicitly. We order the m pixels of the density map 
and form a data vector, x = [(5(no), (5(ni), 5(nm-i)], and 
write the covariance of the data as dj = (xiXj). The likeli- 
hood function is, 

1 



^(27r)'"detC 



exp 



(4) 



The covariance between the overdensity in pixels i and j 
separated by an angle 6ij is given by the sum of the signal 
and the noise components. 



2/4-1 

Cij = -^Viicos e,,)BfCi + N,j 



(5) 



where Nij is the noise covariance matrix and Vi are Legendre 
polynomials. The noise matrix is taken to be diagonal with 
Poisson elements given by Na — wl/n. The finite resolu- 
tion of the pixelised map attenuates the power spectrum by 
the pixel window function, Bi, which depends on the pixel 
geometry (Gorski et al. 2005). 

We now derive a power spectrum estimator that max- 
imises the likelihood function, L. The quadratic form was 
introduced by Tegmark (1997) and explicit derivations have 
been given in Dodelson (2003, Ch. 11), Dahlen & Simons 
(2008) and Bond, Jaffe & Knox (1998). We give an overview 
here, since many variations exist. 

We denote the set of parameters to be estimated by the 
vector A. For our study, \i will represent a bin of the power 
spectrum. We begin with an initial estimate, A'-'^' , and intend 
to use an optimisation algorithm to find a better estimate. A, 
that maximises the likelihood function. With the assumption 
that In L has a quadratic form near the peak, we may apply 
the Newton-Raphson root-finding method to move toward 
the peak of the likelihood function (Press et al. 1992), with. 



A = A 



(0) dlnL/d\ 



lnL/9A2 



(6) 



A(0) 



This expression may be used iteratively to locate the peak. 

We evaluate the derivative terms in Appendix A, and 
now simply state the final result for one iteration step: 



A. 



,N)} 



(7) 



(8) 



The matrix A is a mixing matrix that sets the normalisation 
and may be specified to form linear combinations of the 
bin estimates. We will use this matrix to shape the window 
functions. The second term in Eq. 7 subtracts the noise bias. 

We see that the estimator weights the data by their 
covariance: C~^x. This approach has the favourable prop- 
erty that spatial modes contribute to the measurement with 
an inverse- variance weight. The weighting also appropriately 
'tapers' the map near the mask boundary giving compact 
window functions in harmonic space. 

We have not yet specified the parameter vector A. We 
set A to bins of the three-dimensional power spectrum and 
evaluate the derivative matrix in Eq. 8, as. 



dC, 



dXk dPk 



a — ST 



21 + 1 
An 



Vi (cos 9ij)Bfgi (fc) Ahi k ■ 



(9) 



Here, we have replaced the integral in Limber's equation 
(Eq. 2) with a discrete sum over Infc with logarithmic bin 
width Ain k ■ 

The expectation of the estimate is given by. 



(A) = AFA, 

where we have introduced the Fisher matrix, 



(10) 
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The variance of the estimate is, 



Var(A, A) = AFA' 



(11) 



(12) 



and the window functions are W = AF. 

The inverse of the Fisher matrix represents the min- 
imum variance that we may hope to achieve on A. With 
A — F^^ we see that we have an estimator that is optimal 
in the sense that it is unbiased and has minimum variance 
(Tegmark 1997). This approach may not be practical how- 
ever, because the Fisher matrix is often singular or numeri- 
cally ill-conditioned. Intuitively, this reflects the fundamen- 
tal limit that we cannot probe the power spectrum at scales 
smaller than ~ (AS)^^ where AS is the angular size of 
the survey. 

Instead, we choose A with the aim of diagonalising the 
covariance matrix. By factoring the Fisher matrix as F = 
MM"^, we can set A = M~^. The covariance matrix is 
now Var(A, A) — M~^FM~^ . In practice, we compute M 
as the square root of the Fisher matrix using a singular 
value decomposition (SVD) method. We also rescale M to 
normalise the window functions such that J]]^ Wij = 1. This 
approach was shown by Tegmark (1997) to result in sharper 
window functions than the common choice for A, a diagonal 

matrix with An = Fij^ 

We find that the matrix M is ill-conditioned when the 
window functions are broad, especially for the SV sample 
which has a wide redshift distribution. To find a stable inver- 
sion, we use a pseudo-inverse technique by keeping only the 
largest singular values. The consequence of using a pseudo 
inverse is that the covariance matrix will not be perfectly di- 
agonalised. The covariance matrix for the VIPERS estimate 
after carrying out this operation is shown in Fig. 3. The 
choice of how many modes to keep in the pseudo-inverse af- 
fects the shape of the window functions and the scales that 
are probed. We find that the smaller singular values probe 
large scales, in the same fashion as in the SVD analyses by 
Eisenstein & Zaldarriaga (2001) and Mailer et al. (2005). 
We set the scale ensuring that the resulting window func- 
tions are positive and reach the largest scales available to 
the survey. 

The estimate and covariance model depend on the cho- 
sen fiducial cosmology through the form of the likelihood 
function and the projection kernel. Both the shape and nor- 
malisation of the power spectrum can be important. Al- 
though the normalisation cancels in the estimator (neglect- 
ing the noise term), it is important for the Fisher matrix 
and covariance. For these reasons, maximum likelihood es- 
timators are often applied iteratively to arrive at consistent 
results. We explore these dependencies with simulations in 
the next section. 



5 SIMULATIONS 

5.1 Gaussian realisations 

As a test of the method, we estimate the power spectrum 
for Gaussian realisations of the projected density field. The 
simulations are constructed using a CAMB power spectrum 
with the Halofit model (Lewis, Challinor & Lasenby 2000; 
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Figure 4. The top frame shows the recovered power spectrum 
from the mean of 1000 independent Gaussian simulations. The 
theory is convolved with the window functions (plotted at bot- 
tom) and we find that it matches the measurement to within a 
few percent. The error bars computed analytically from the Fisher 
matrix (shaded area) agree with the distribution of MC runs (out- 
line). The noise term in Eq. 7 is shown as a dashed curve. 
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Figure 5. In the top panel we plot the projection kernels, gi{k). 
We compare the kernels derived from the Wl redshift distribu- 
tion only (solid curves) and those from W4 only (dashed curves) . 
Lower panel: the window functions found for the 21deg^ simula- 
tion field are plotted. To check robustness we again compare the 
results derived from the Wl and W4 fields separately. The sec- 
ond peak in the window functions at A; > 1 h Mpc~^ arises from 
the pixel scale of the map; beyond fc = 2 h Mpc~^ the window 
functions rapidly drop to 0. 
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Figure 6. The influence of tiie assumed redsiiift distribution on 
the deprojection. Top panels: the simulations are analysed with 
the Wl and W4 redshift distributions. Bottom panels: for an 
overly-conservative test, we modify the VIPERS redshift distri- 
bution to adjust the mean redshift by Az = ±0.05 (labelled Test 
1 and Test 2) leading to systematic shifts in the amplitude of the 
estimated power spectrum of -|-6% and —10%. The redshift dis- 
tributions are plotted in the right panels and the derived power 
spectra are on the left. In the top and bottom, the filled grey dis- 
tribution represents the redshift distribution of the full VIPERS 
sample. 



Smith et al. 2003) and projected through the VIPERS red- 
shift distribution. With these Ci , we use Healpix Synfast 
to generate simulated density maps. We produced 1000 in- 
dependent maps with a resolution of 7' (nsidc = 512). We 
add noise fluctuations by drawing from a Gaussian distri- 
bution assuming a variance given by Poisson statistics with 
n = 50 galaxies/cell. This is a higher level of noise than 
we find in the CFHTLS catalogue. For the geometry of the 
mock survey, we use the actual survey mask of the W2 field. 
It includes 1592 pixels covering 21 deg^. 

We compute the Fisher matrix in logarithmic bins from 
fc = 0.01 - 100 with Aiogfc = 0.05. The resuhs are plot- 
ted in Fig. 4. The data points are re-binned to Aiogfc =0.1 
and plotted from k = 0.06 — 0.7 h Mpc~^. The correspond- 
ing window functions are shown in the bottom panel. The 
data are plotted at the peaks of the window functions. On 
large and small scales, the window functions begin to over- 
lap and converge as the limits set by the survey geometry 
are reached. On small scales, we see a secondary peak in the 
window function at ~ 2 h Mpc~^ which arises from the 
pixel scale of the map, see Fig. 5. 

We convolve the theory power spectrum with the win- 
dow functions and find that the mean of the Monte Carlo 
runs agrees well within a few percent. We expect that the 
precision is limited by the finite binning of the Fisher matrix 
and truncation of the window functions, but these effects are 
well below the statistical uncertainties. The errors computed 
analytically from the Fisher matrix agree with the distribu- 
tion of Monte Carlo runs to within a few percent. These 
errors are for a single field, and so, we can expect to achieve 
a factor of two better with the combination of four fields. 
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Figure 7. The sensitivity of the estimator to the assumed fidu- 
cial model, (a) We vary f2m keeping other parameters fixed. The 
top frame shows the fiducial power spectra (solid lines) and the 
points mark the derived power spectrum measurements. The bot- 
tom frame shows the percent differences from the reference model 
for each trial. The correct shape is recovered, but there is a shift 
in the amplitude of the estimate due to the geometric dependence 
of the projection kernel on the cosmology, (b) We modulate the 
small-scale amplitude of the fiducial power spectrum with an in- 
terpolation parameter a„; ; a„; = and 1 correspond to the linear 
and Halofit models, respectively. We find that the derived power 
spectrum is not sensitive to the small-scale amplitude. 



5.2 Dependence on redshift distribution 

We checked the robustness of the measurement to uncer- 
tainties in the redshift distribution by repeating the analysis 
with different assumed distributions. The simulations were 
generated with the measured redshift distribution from the 
complete VIPERS sample, and we first re-analysed them 
with distributions derived from two subsamples, the redshift 
distribution of Wl and W4. The numbers of spectra taken in 
the two fields are similar and the spectroscopy covers similar 
areas, but, the fields are widely separated on the sky and so 
any differences could be attributed to cosmic variance. We 
compute the Fisher matrix using the two distributions and 
find that the projection kernels and window functions agree 
(Fig. 5). The bias introduced by a mismatched redshift dis- 
tribution is at the percent level, below the statistical errors. 
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Additionally, the measured redshift distribution could 
be inaccurate due to sampling biases in the VIPERS sur- 
vey. In Section 2.3, we concluded that the uncertainty in 
the mean redshift of the distribution, z, is known to better 
than Az = 0.01. As an overly conservative check, we exam- 
ine the consequences of shifting the redshift distribution by 
Az = ±0.05. This was done by modulating the measured 
distribution of the full VIPERS sample by the linear func- 
tion /(z) = 1 ± 1.5(z — z). The modified distributions have 
zi = 0.656 and Z2 = 0.752, while the original sample has 
z = 0.703, see the lower panel of Fig. 6. We find that reduc- 
ing z by 7% lowered the derived power spectrum by 10%. 
Increasing z by 7% increased the power spectrum by 6%. 
This large shift in z would thus lead to a systematic error in 
the estimated bias factor at the level of 3 — 5%. 

5.3 Dependence on fiducial cosmology 

The dependence on the fiducial cosmology enters the anal- 
ysis in two ways. First, we rely on the cosmology to model 
the likelihood function. In the maximum likelihood estima- 
tor, the data covariancc matrix plays the role of a weight. 
Modifying the fiducial power spectrum changes the weight- 
ing function and could bias the estimator. We can expect 
that assuming the wrong matter density for example, could 
bias the estimator and make the variance properties sub- 
optimal. 

The second dependence on the fiducial cosmology is 
through the projection kernel. In the previous section, we 
discussed how shifting the redshift distribution affects the 
amplitude of the power spectrum estimate. We can expect 
that varying the cosmology and the redshift-distance rela- 
tion will have a similar effect. 

Our Gaussian simulations were constructed using the 
reference ACDM power spectrum with the Halofit model. 
To test the dependence on the cosmology, we first re-analyse 
the maps using fiducial power spectra with different assumed 
values of the matter density, taking Q,m = 0.25, 0.30 and 
0.35. All other parameters were held fixed at the reference 
values. We find that despite assuming the wrong matter den- 
sity, we recover the correct shape of the power spectrum from 
the mean of 1000 simulation runs to within 2%, see Fig. 7, 
Panel a. However, it is clear that the amplitude is strongly 
biased. This is due to the dependence of the projection ker- 
nels on Q,m- This geometric dependence on the background 
cosmology dominates over any bias in the estimator due to 
sub-optimal weighting. These findings support an iterative 
approach. 

Next, we check the influence of variations in the ampli- 
tude of the power spectrum at small scales. The shape of 
the power spectrum on small scales has developed with the 
aid of n-body simulations but it remains a source of system- 
atic uncertainty. We vary the small-scale amplitude using an 
interpolation parameter a„;: 

Pk = Pk,Hn + {Pk,nl — Pk,lin) dnl ■ (13) 

We test a range of amplitudes with a-ni = {0,0.5,1,2} 
(a„; = 1 gives the Halofit model). We find that the estimator 
is remarkably robust, see Fig. 7, Panel b. The discrepancy 
introduced by the variation in the small-scale amplitude is 
less than 2% on large scales and it is dominated by numeri- 
cal uncertainties up to ~ 0.2 h Mpc~^. This supports the 



conclusion that using sub-optimal weights does not signifi- 
cantly bias the result. 



6 RESULTS 

We now carry out the estimation of the power spectrum on 
the CFHTLS data for the VlPERS-like sample and three 
photometric redshift subsamples, see Fig. 8. For the analy- 
sis, a fiducial CAMB power spectrum with the Halofit model 
is assumed. The Fisher matrix is computed with 60 bins log- 
arithmically spaced from k — 0.01 — 10 with Alogfc = 0.05. 
We use a wide fc-range to map out the window functions 
but all these data are not useful for analysis. We restrict the 
study to 13 points from k = 0.03 — 0.6 h Mpc"'^. We can go 
to smaller scales, although the Gaussian error estimate will 
not be appropriate. We use a bin size of A log fe = 0.1 which 
is appropriate choice considering the width of the survey 
window functions. The plotted error bars are derived from 
the diagonal elements of the covariance matrix found com- 
puted from Eq. 12. 

For each field, we compute the normalised quantity from 
Eq. 7, Dj = I {x^EjX-Tr(EjN)}, where k indexes the 
fields 1-4, along with the Fisher matrix (Eq. 11). These re- 
sults are then summed together, and the final combined 
estimate is computed by Ai = | 'I2t=i 5Zj ^'iVj where 

A = {J2t=i F'')'^^'^ normalised such that J2i i^F)^. = 1. 
This combination properly weights the data. The covariance 
of the estimate for the VIPERS-like sample is shown in Fig. 
3. At low k, neighbouring bins are nearly 50% correlated, 
but the matrix becomes more diagonal at larger k. The lim- 
itation in diagonalising the covariance matrix comes in the 
inversion of M = F^^^. This is computed with a pseudo- 
inverse method. The inversion becomes easier for the nar- 
rower photometric redshift slices where a nearly perfect in- 
version is possible. The window functions are sharper for 
these redshift slices as well. 

We do not run the maximum likelihood algorithm in an 
iterative fashion. The data do not support strong constraints 
on ACDM parameters alone and we find that beginning with 
a fiducial ACDM power spectrum, we have a very good fit 
to the data. This indicates that our starting point is already 
near to the peak of the (very broad) likelihood function. 
However, we do effectively carry out one iteration of the 
estimator to find the galaxy bias and set the amplitude of 
the fiducial power spectrum. This is necessary because the 
estimator and covariance do not simply scale linearly with 
amplitude in the presence of noise. A second run allows us 
to set the amplitude of the fiducial power spectrum ensuring 
that the error estimate is correct. 

We compute a one-parameter fit to estimate the galaxy 
bias on linear scales. We restrict this fit to the first 9 points 
at fc < 0.2 h Mpc~^. Given the Pk measurements in vector 
d and the convolved ACDM model in vector m, we find the 
amplitude, o, that maximises the likelihood function. In L = 
—1/2 (d — am)^ ^kk' ("^ ~ am). The solution is given by, 



with variance al = ( m.^ C^^/Hij . The resulting values of 
the galaxy bias are listed in Table 3, where we have assumed 
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Figure 8. The deprojected three-dimensional power spectrum in logarithmic bins for the VIPERS-like sample and three photometric 
redshift selected samples: Zpi^^t = 0.5 — 0.6, 0.6 — 0.8 and 0.8 — 1.0. The window function of each band is shown in the panel below each 
plot. The theoretical power spectrum (both linear and Halofit models) is convolved with the window function and overplotted with the 
best-fit linear bias listed in Table 3. We use the first 9 data points, up to fc ~ 0.2 h Mpc~^ indicated by the vertical dotted line, to 
estimate the linear bias. The corresponding comoving distance and angular scale (at z = 0.7) are included as a guide. 



a value of — 0.8. The bias increases with redshift as ex- 
pected for a flux limited survey. In fact, the amplitude of 
the power spectrum is not seen to change with increasing 
redshift, indicating that the evolution of the growth factor 
and the galaxy bias factors approximately cancel. We find 
an approximately constant error on the bias factor in each 
redshift range; this is simply due to the fact that the ampli- 
tude of the fiducial power spectrum (from which the errors 
are derived) is approximately constant. 

Also in Table 3, we list the values at the best fit. The 
number of degrees of freedom is approximately 8. The 
values are lower than expected, specifically for S8 for which 
we find = 2. Formally, the probability of finding ^ 2 
with 8 degrees of freedom is 0.019. This could indicate that 
the covariances and, consequently the error bars, are over- 
estimated for this sample. 



7 PARAMETER ESTIMATION 
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Figure 9. The covariance between the photometric redshift sam- 
ples is plotted. The circle markers show the full computation of 
the variance of the S7 slice that accounts for the survey geometry. 
The solid curves show the analytic Fisher approximation in the 
full-sky limit (Eq. 17). Inset is the full correlation matrix for the 
three samples. 



Before we may carry out a joint analysis, we must estimate 
the covariance between the overlapping photometric sam- 
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Table 3. Best-fit galaxy bias 



Sample 


69 




VP: VIPERS-like 


1.38 ±0.05 


7.2 


S6: 0.5 < Zphot < 0.6 


1.36 ±0.05 


7.1 


S7: 0.6 < Zphot < 0.8 


1.52 ±0.05 


5.7 


S8: 0.8 < Zphot < 1-0 


1.68 ±0.05 


2.0 



Table 4. Marginalised parameter estimates 





0.30 ±0.06 


bs6 


1.39 ±0.08 


bsT 


1.55 ±0.08 


bss 


1.72 ±0.10 



pies. We will estimate the covariance of the estimators with a 
Fisher matrix approach in the full-sky limit and then rescale 
to find the errors for our survey geometry. 

We compute the covariance between two samples la- 
belled A and B. To simplify the expression, we write the 
product of the kernel with the beam and integration step 
as, gi{k) = gi{k)Bf Aink and the sum of the signal and noise 
covariance as, C; = CiBf + 4p. It is more convenient to use 
the harmonic space representation, and we switch the data 
vector from Xi to aim- The covariance matrix is diagonal: 

ClmWrn' = SuiSmm'Cl- 

A component of the Fisher matrix for a single sample 
A for two power spectrum bins, k and k' is. 

Fa,,,' = E ^-^-mikygi'ik') (cr) (15) 

I 

We can write the quadratic estimator for the power spec- 
trum as. 



(16) 



We find the covariance between the two sample esti- 
mates to be, 



-mU)9f{k') 



(17) 



We scale by the fractional sky coverage of the survey, faky, 
which approximately accounts for the number of modes that 
may be probed. The small survey size also broadens the 
window functions, which we account for in the covariance 
with: C',,, = W^Cfefc,WS^. 

The variances for one sample computed with Eq. 17 in 
the full-sky limit match well with the full computation of the 
Fisher matrix (Eq. 12), see Fig. 9. Although, we find that 
the full-sky computation underestimates the variance by a 
factor of ~ 2. This is not surprising since we have neglected 
the precise survey geometry. As a correction, we rescale the 
estimate to match the variance in the S7 slice. In Fig. 9 we 
also show the analytic estimates of the covariances between 
the three redshift slices that we may now use to perform a 
joint likelihood analysis. 

Using the sample covariances, we jointly estimate the 
linear galaxy bias factors of the three photometric redshift 
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Figure 10. The joint likelihood surfaces of Q,m and the bias 
parameters for the three photo-z samples (6s6i bsT, bsg). The 
inner and outer contours indicate the 68% and 95% confidence 
level. The marginalised likelihoods of each parameter are listed 
in Table 4. The fit is limited to the first 9 data points, giving 
kmax = 0.2 h Mpc-i. 



slices, S6, S7 and S8, labelled as bse, bs7 and bss along 
with Qrn- AH other ACDM parameters are held fixed and 
we set (Tg — 0.8. We compute the likelihood of a model with 
the full covariance matrix. The fit is limited to the first 9 
datapoints of each sample, giving a maximum k of kmax ~ 
0.2 h Mpc~^. We exhaustively evaluate the likelihood over 
the 4-dimensional parameter grid. Views of the likelihood 
surface, marginalised over pairs of parameters are shown in 
Fig. 10. The marginalised constraints are listed in Table 4 
with 68% confidence intervals. 

The joint analysis prefers a slightly higher value of flm , 
0.30 ± 0.06, versus the fiducial model with 0.272. Due to 
the correlations between parameters, this results in higher 
values of the galaxy bias factors than were found with the 
fiducial model fixed (Table 3). 



8 CONCLUSIONS 

The CFHTLS-Wide fields probe a significant cosmological 
volume at redshifts not reached by other galaxy surveys to 
date. We use the projected density field from photomet- 
ric redshift samples to constrain the real-space power spec- 
trum and derive constraints on the matter density and lin- 
ear galaxy bias factors. These results are made possible by 
precise knowledge of the redshift distributions provided by 
preliminary results from the VIPERS survey. 

The primary advantage of computing the power spec- 
trum directly from the angular distribution, instead of us- 
ing conventional spherical harmonics Ci is that we may con- 
struct window functions in Fourier space. By optimising this, 
we achieve sharper constraints on the power spectrum than 
when we are limited to £ bands. This approach comes with 
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the cost that we must adopt a fiducial power spectrum. We 
showed that using the wrong fiducial power spectrum, al- 
though leading to sub-optimal weights, does not significantly 
bias the estimate. This is true even on small scales, and we 
can effectively use this method to deconvolve small and large 
scales in Limber's equation. Residual systematic error on the 
derived power spectrum is at the 1% level, well below the 
sensitivity of the measurement. 

The deprojection does strongly depend on the assumed 
redshift distribution of the galaxy sample as well as the cos- 
mology used to compute the redshift-distance relation. The 
cosmology dependence of the measurement makes the in- 
terpretation difficult, but to a first approximation, only the 
amplitude is affected; the shape of the power spectrum is 
recovered correctly. Thus, a converging iterative procedure 
can be implemented by updating the fiducial model and re- 
peating the analysis. 

There is a degeneracy between a shift in the assumed 
redshift distribution and the cosmological model. This is un- 
avoidable when studying a field in projection. However, the 
constraints on the redshift distribution can always be im- 
proved with further observation. In our analysis, from the 
sampling biases present in the VIPERS spectroscopy, we es- 
timate the uncertainty in the mean redshift to be at the 
1% level. Thus, we do not expect a strong systematic er- 
ror in the derived galaxy bias parameters. We do note that 
the observed trend of low values for the best-fit models 
in the higher redshift samples can arise if the covariance is 
overestimated. This could be a weak hint that the true mean 
redshift is lower than what we assume or that a modification 
is needed in the fiducial cosmology. 

Recently, the galaxy bias was measured from the 
CFHTLS-Wide fields in the context of the halo model by 
Coupon et al. (2011). Our final two photometric redshift 
bins, S7 and S8, correspond with samples constructed by 
Coupon et al. so we are able to compare the resulting bias 
values. Coupon et al. constructed volume-limited samples 
using luminosity cuts resulting in a selection of brighter 
galaxies, thus we may expect their bias values to be larger. 
The halo model constraints of Coupon et al. (2011) give for 
S7, bg = 1.44 ±0.01, and for S8, bg = 1.79 ±0.03. These val- 
ues have been scaled by 1.03 to transform from a cosmology 
with Qrn = 0.25 to flm ~ 0.272 which is assumed here. Our 
value of bg for the S7 sample is higher, while for the S8 sam- 
ple it is lower, although both are in agreement with Coupon 
et al. within the 2cr confidence limit. The measurements are 
based on different physical scales (Coupon et al. restrict the 
correlation function to angular scales < 1.5°) and different 
model assumptions have been used. Thus it is reasonable to 
consider the measurements as independent estimates. 

Our results provide a preliminary look at the large-scale 
structure field probed by the VIPERS colour selection and 
demonstrate the strengths of the VIPERS sample for clus- 
tering studies at z > 0.5. We anticipate promising results 
with the full VIPERS spectroscopic sample. 
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APPENDIX A: QUADRATIC ESTIMATOR 



In Section 4, we reasoned that we would apply the Newton-Raphson root-finding algorithm (Eq. 6) to locate the peak of the 
likelihood function (Eq. 4). We now continue and evaluate the derivatives of the likelihood function, ^^^^ and 

assume that the covariance of the data depends linearly on the parameters: Cij = Pk,ij\k + Nij. 
The first derivative term is, 



9' '"^ We 



ainL 



aindetC t9C-i 

+ X ^7:: — X 



dXi 



dXi 



Two identities have been used: In(detC) = Tr(lnC) and = — C~^^C~^. The second derivative, or curvature, is 



dXidXi \ dXi dX-i 



+ 2x^C-i|^C-i|^C-ix 

oXi aXj 



(Al) 
(A2) 

(A3) 



We neglect the second derivative terms. To simplify, we replace the curvature by its average over an ensemble of realizations 
of the data. This is known as the Fisher matrix. 



1 InL 

2^dX,dXj' 



-Trie 



dXi dXj 



(A4) 



Inserting these expressions into Eq. 6, we find that one iteration step in the Newton-Raphson algorithm is given by. 



AO) 



Z^{dXidxJ 



dlnL 
dXi 



d'^ InL 



dXidXi 



a: 



(0) 



x^C 



dXi 



C"^x-Tr ( C 



_i9C 
dXi 



(A5) 
(A6) 

(A7) 



The terms on the right are computed with the parameter set A^°' . We may simplify further by rewriting the trace term with. 



dXi 



oXi 



(A8) 

(A9) 

(AlO) 



where we use the linear dependence of C on the parameters. Substituting into A7, the product of the Fisher matrix with its 
inverse leads to a cancellation of the A^°' terms. We are left with the final estimator in quadratic form. 



Ai - 9 XI ^ 
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■ i^C 
dXi 
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dXi 



(All) 
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